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While not generally a conservation law, any symmetry of the equations of motion implies a 
useful reduction of any second-order equationto a first-order equation between invariants, whose 
solutions (first integrals) can then be integrated by quadrature (Lie's Theorem on the solvability 
of differential equations). We illustrate this theorem by applying scale invariance to the equations 
for the hydrostatic equilibrium of stars in local thermodynamic equilibrium: Scaling symmetry 
reduces the Lane-Emden equation to a first-order equation between scale invariants u n ,v n , whose 
phase diagram encapsulate all the properties of index-n polytropes. From this reduced equation, we 
obtain the regular (Emden) solutions and demonstrate graphically how they transform under scale 
O . transformations. 
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I. SYMMETRY REDUCES THE ORDER OF ANY SECOND ORDER ODE 

Lie showed how the invariance of a second-order ordinary differential equation (ODE) under a point symmetry leads 
constructively to a first-order ODE plus a quadrature. If the symmetry were a variational symmetry of the Action 
5_i 1 Principle, Noether's well-known theorem would lead to a conservation law [H-Q- 

We consider the scaling symmetry £ — > At;, #„(£) — » 9 n A(0 of the second-order ODE Lane-Emden equation (Section 
III) which describes the hydrostatic equilibrium of a gaseous sphere or star in local thermodynamic equilibrium (Section 
II) . Because this is not a variational symmetry, but only a symmetry of the Lane-Emden equation, scaling symmetry 
leads only to a non-conservation law [4|, [5|. This is a first-order differential equation for scale invariants u n , v n , 
which can be solved for v n (u n ) for given boundary conditions (Section III). From these first integrals, quadrature 
finally leads to solutions of the original second-order equation (Section IV). 

Sections II and III will consider only regular solutions which have density finite at the origin and apply only 
to complete polytropes. Section IV will generalize to irregular solutions which have densities infinite at the origin 
(F-solutions) or vanish away from the origin (M-solutions) and can apply only to stellar envelopes. 
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FIG. 1: Emden Functions (Regular Solutions of the Lane-Emden Equation) of order n=l,2,3,4,5 

II. HYDROSTATIC EQUILIBRIUM OF SELF-GRAVITATING SPHERES 

An adiabatic sphere in hydrostatic equilibrium obeys the equations of equilibrium between gravitational and internal 
(pressure)forces and of mass continuity 

— dP/pdr = Gm/r 2 , dm/dr = 47rr 2 p , (1) 

where the local pressure, mass density, and included mass P(r), p(r), m(r) depend on radius r. In terms of the the 
gravitational potential V(r) = J Gm/r 2 dr and thermostatic potential (specific enthalpy, ejection energy) H(r) = 

Jq dP/p, (P) and its integrated form 

CM 

- dH/dr = dV/dr, V(r) + H(r) = — , (2) 

it 

expresses the conservation of the specific energy as sum of gravitational and internal energies, in a star of mass M 
and radius R. The two first-order equations ((T|) can be combined into the second-order Poisson's Law 

-Mr 2d -^-)+^Gp{H) = Q , (3) 

r z dr \ dr I 

in terms of the enthalpy H{r). 

The two equations ([1]) can always be written 

d\og u/dlogr = 3 — u(r) — n(r)v(r) , d log w/d log r = u — 1 + v(r) — d log [1 + n(r)]/d log r , (4) 

in terms of the logarithmic derivatives 

u{r) :— d log m/d log rw(r) := n(r)v(r) — — d log p/d log rv(r) :— — d log (P/p) /d log r (5) 

and an index n(r) 

n(r) := d log p/d log (P/p) , 1 + l/n(r) := d log P/d log p , (6) 
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FIG. 2: 'Mass' included inside radius £ for complete polytropes of order n=l,2,3,4,5 terminating at fi n = 
3.141,4.353,6.897,14.97,oo. 

which depends on the local adiabatic equation of state P — P(p). 

For solutions regular at the origin, spherical symmetry requires that dP/dr — and that p{r), P(r), H{r) be even 
functions of r. Mass continuity requires, to order r 2 , 

p{r)* Pc {l-Ar% m{r)^^.{l-\A^)^^.pTp{rf' 5 . (7) 

3 5 3 

Thus, near the origin, the average mass density inside radius p(r) :— 4 ™Jjj 3 w pc^ 5 p{r) 3 ^ 5 and u(0) = 3,u(0) = 0. 
At the stellar surface, u(R) = 0, v(R) = oo. Here p(r) := 4^3/3 is the average mass density inside radius, so 
that u(r) = 3p/p decreases from u(0) = 3 at the origin to u(R) = oo at the stellar surface. The ratio v(r) := 
|(— Gm/r)/(P/p) = | ('gravitational energy')/('internal energy of ideal gas') increases from v(0) = at the origin to 
v(R) = oo at the stellar surface. 

III. SCALE TRANSFORMATIONS ON SOLUTIONS OF THE LANE-EMDEN EQUATION 

The symmetry we consider is scaling symmetry, the most general simplification that one can make for any dynamical 
system. If a scale transformation r — > Ar transforms m, p, P multiplicatively, the logarithmic derivatives u, v will 
be scale invariant. The structural equations (5) will then be autonomous, if and only if n = constant, so that, 
P(r) = Kp(r) 1+ ™ , with the constant K determined by the constant specific entropy. When this is so, scaling 
symmetry leads to the two coupled autonomous equations for the scale invariants 

du n /dlogr — u„(3 - u n - nv n ) , dv n /dlogr = v n {u n — 1+ v n ) , (8) 

implying the first-order equation 

dv n /du n = v n (u n - 1 +v n )/u n (3 —u n - nv n ) . (9) 

This reduced equation for v n {u n ) incorporates all the consequences of scaling symmetry, from which all the properties 
of polytropes follow. 
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FIG. 3: Emden function #3(£, A) for three scales A=4 (blue), 1 (green), .25 (red). On a log-log plot, rescaling simply transforms 
the Lane-Emden functions along lines of constant 113 := —d log (P/p)/d\ogr. The three dashed lines connect homologous points 
for the three illustrative values V3 = 0.081, 0.25, 2.3. 



In this and the next section, we consider only regular solutions (Emden or E-solutions) of the Lane-Emden equation, 
which have finite density at the origin, and are applicable only to complete polytropes. At the origin, the initial 
conditions on regular solutions of [9] are u„(0) = 3, i> n (0) = 0. (We defer to Section IV the irregular solutions v n {u n ), 
where the density at the origin is not finite.) 

Figure 3 shows how, moving radially outwards, u n (r) = ip{r) / p{r) decreases from 3 at the origin to at the outer 
boundary R and v n (r) increases from to 00. The combinations u) n {r) := (to") 1 """ 1 ' = — £, 1+Un 9' n approach the 
finite values oui n (R) and characterize each n-polytrope. 

In terms of the dimensional central density, pressure p c , P c = Kpl +1 ^ n and constant a 2 := ((71 + 
l)/4ftG)Kpl^ n ~ 1 , H c := (n + l)(P/p) c = (n + l)Kp X J n , the dimensional radius, enthalpy, mass density and in- 
cluded mass are 

r = a£, H = H c 6 ni p = p c 6%, m(r) - (4 W 3 )(-^X) . (10) 
where prime designates the derivative ' := d/d^. The scale invariants are 

u n :=-ie n n /6' n , v n x=-£#J9 n u n := {uv n fl^ = , (11) 

where uJ n := 2/ (n — 1). 

In these dimensionless units, , Poisson's Law ([3]) becomes the Lane-Emden equation Q 



(12) 



whose regular solutions have 

9 n = const, 0' n (0) = 0, u(0) = 3, v n (0) = , (13) 

The normalized regular solutions with 9 n = 1 define the Emden Functions for the n-polytrope radial evolution (Figure 
1). Their dimensionless normalized 'mass' —0 6' n (^) is shown in Figure 2. The first order equations (jSJ have analytic 
solutions listed in Table 1 for n = 0, 1, 5, but must be integrated numerically for any other polytropic index n. 
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FIG. 4: 'Mass' included inside radius £ for n=3 polytrope. On a log-log plot, rescaling simply transforms the 'mass' along 
lines of constant V3 := — d log [P/p)/d log r. The dashed lines connect homologous points for the three illustrative values 
v 3 = 0.081,0.25,2.3. 



The Lane-Emden equation (TT21 is invariant under the scale transformation £ — )• ^4£, t — > t + log A, — >• 
A u ' n 8 n (A £) := 6 u a{£)- Besides the Emden Function 9 n (^), normalized so that 9 n i(£) = 0n(Q, the Lane-Emden 
equation has rescaled regular solutions # nJ 4(£), whose value at the origin is A Un . Figures 3 and 4 show log- log plots 
of the rescaled n=3 function 03a(O and rescaled 'mass' —9' 2, &^9za{0 for three different rescalings A=4 (blue), 1 
(green), 0.25 (red). On a log-log plot, rescaling appears as a translation along the dashed lines of constant v n and 
u n , which connect homologous points. All the familiar properties of polytropes Q follow from this scaling symmetry 
incorporated in the reduced equation for v n (u n ). 



TABLE I: Scaling Invariants and normalized Emden Functions for n=0, 1, 5 
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V n (u n ) 
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"n(0 


0»(fl 


„ , (n+i) 


fin 


0^n 





u — 3 
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2*7(6 -O 


l-C 2 /6 
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0.333 
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parametric u„(f ),««(£) 


^/(l-^cotf) 
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3.14 
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1 - u 5 /3 
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£7(3 + £ 2 ) 


(l + C 2 /3)- 1/2 
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00 






IV. PHASE DIAGRAM FOR ALL FIRST INTEGRALS 



Equation © and the integral curves v n (u n ) and u n (v n ) are scale invariant functions of the scale invariant (r/R) 
plotted in Figure 5 for n=2, 3, 4, 4.99. All the dependence on scale R is contained in ©which reads 

Jo V n (U n (V n ) - 1 + Vn) J 3 U n (S - U n - nVn(U n )) 

for the regular solutions obeying u n (0) = 3, v n (0) — at the origin. 

The irregular solutions of © apply only to incomplete polytropes, the scale-invariant envelopes of stars. As shown 
in Figure 6, these are F-solutions (green curves, u n (0) > 3) which have infinite density at the origin and M-solutions 
(blue curves, v n (fi)) which have density vanishing before the origin (v(0) > 0). These irregular solutions would obtain 
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FIG. 5: The Emden scale invariants u„(r) = dlog m/dlog r, v n (r) = — d log P/dlog r and u) n (r) = (uv n ) 1 ^"~ 1 K While v„(r) 
diverges at the stellar radius R, u> n (in green) approaches the finite values qoj2 = 10.49, o^3 = 2.02, o^4 = 0.73, 0^4.99 ~ 
characteristic of each n-polytrope. The curves in the n = 4.99 figure are excellent approximations to the n = 5 functions on 
the bottom line of Table I for which R = 00. 

by integrating in from the boundary values F^n >o w n and ouj n >m w„. Their separatrices (red curves) are the 
regular (Emden) solutions, which have finite density at the origin, u n (0) — 3, v n (0) = and boundary values ow n . 

V. CONCLUDING SUMMARY 

A symmetry of the equations of motion generally implies, not a conservation law, but a still useful reduction of order 
of the equations of motion to a first-order equation between invariants, which can then be integrated by quadrature. 

For the Lane-Emden equation, the reduced equation is the first-order equation ^ , whose first integrals (Figure 6) 
encapsulate all the properties of index-n polytropes. From this reduced equation (Figure 5), we obtained the regular 
(Emden) solutions (Figures 1, 2) and simply demonstrated their scale dependence on log- log plots (Figures 3, 4). 
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FIG. 6: Critical Points and Phase Diagrams for v n (un) for n=2, 3, 4, 5. The green curves (F-solutions) have infinite density at 
the stellar center. The blue curves (M-solutions) have vanishing density away from the origin. Their separatrix, the red curve 
(Emden solutions), has finite density at the origin. 
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